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Abstract 

We discuss the lightest scalar resonances, /o(600), k(800), ao(980) and 
/o(980) in the extended Nambu Jona-Lasinio model. We find that the model 
parameters can be tuned, but unnaturally, to accommodate those scalars ex- 
cept the /o(980). We also discuss problems encountered in the K Matrix 
unitarization approximation by using N c counting technique. 
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1 Introduction 

The original model of Nambu and Jona-Lasinio [1] (NJL)was proposed as a dynami- 
cal model of the strong interactions between nucleons and pions, before the invention 
of QCD. In this model pions appear as the massless composite bosons associated 
with the dynamical spontaneous breakdown of the chiral symmetry of the initial 
lagrangian. Even after the invention of QCD, the NJL model or the extended NJL 
(ENJL) model still serves as a useful tool widely discussed in the literature when 
discussing low energy strong interaction physics at the quark-level, starting from 
last seventies and eighties [2, 3]. More recent extensive reviews can be found in [4]. 

The ENJL model provides a natural extension and hence is considered more 
general than the linear sigma model or the model in which the p meson takes the 
role of a massive gauge boson of the isospin symmetry. The ENJL model attempts to 
provide a unified description to both the scalar sector as well as the vector sector, in a 
chiral symmetric way. However, there have been controversies for a long time on the 
spectrum of the lowest lying scalar nonet in strong interactions. This situation is also 
reflected in the early studies on ENJL models. In Ref. [5], the ENJL model is used to 
study the lightest scalar nonet, and it is found that, m a = 500MeV, m^, = 840MeV, 
fflj ~ 760MeV, m a% = 950MeV. The mass of ao(980) meson (which is degenerate 
to the SU(2) a in Ref. [5]) could not be explained by the NJL model itself, and 
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the mass difference is ascribed to a possible qqqq content in the physical oq. In 
Ref. [6] the role of QCD Ua(1) anomaly is considered and a sum rule is obtained 
between the mass of scalars and the pseudoscalars in the NJL model via 't Hooft's 
instanton interaction [7]. The instanton effect will break the degeneracy between 
the octet and the singlet, which lifts the former and suppresses the latter. Starting 
with a 'bare' quarkonium mass of HOOMeV in Ref. [6], Ua{1) splits the + nonet 
into a singlet ~ lOOOMeV and an octet ~ 1300MeV. Broken SUf(3) further splits 
the masses so that one gets: m ao = 1320MeV, m CTg = 1590MeV, m ao = lOOOMeV, 
rrtK* = 1430MeV. Apparently this assignment is not for the lightest scalars since it 
contains a heavy ao- The above work did not include k (or ^(800)) in their lightest 
scalar nonet. This situation is improved by Volkov et al [8] who discussed the NJL 
model with the 't Hooft interaction (which splits the mass between ao and a) and 
find m ao = 810MeV, m a = 550MeV, m a , = 1130MeV, m K * = 960MeV. Nevertheless 
there is an apparent problem with these results, that is for a small a mass around 
550MeV (and also a small mass for the k), it does not possess enough phase space 
to develop a large width for the a as is revealed by recent determinations. More 
recent work [9, 10] observed that there is growing evidence that a (980), -K"q(800) 
or k, as well as /o(600) or o and /o(980), are members of the low-lying scalar nonet. 
Nevertheless Ref. [9] suffers from the similar problem as in Ref. [8], it also gives a 
rather small mass of a with which it is difficult to explain the large width of the 
sigma simultaneously. 

On the other side, progress has been made in recent few years, demonstrating 
the existence of the light and broad a (or /o(600)) and k resonances [11]. The pole 
locations of a and k are determined using dispersive approaches [12]— [16]. This new 
information on the pole locations urges and enables us to watch more carefully the 
dynamics with respect to the lightest scalars, within the scheme of the extended 
NJL model, which is the purpose of this paper. 

The basic idea of the present paper is outlined already in Ref. [17], where we 
pointed out that in order to understand correctly the mass relations among light- 
est scalars one has to take into account the additional information provided by the 
widths of these scalars, ranging from a few tens MeV to a few hundred MeV. Es- 
pecially when there appears a large width, since it is an unambiguous signal for 
strong interactions in the given channel, the bare mass spectrum at tree level has 
to be strongly distorted. A certain unitarization procedure is necessary to explore 
the relation between the pole mass parameters and the bare mass parameters put 
in the lagrangian. For example it is suggested in Ref. [17] that a a pole locates at 
y/z = m — iT /2 = 470 — 285i MeV corresponds to a bare mass M a ~ 930MeV with 
some uncertainties. 
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This paper is devoted to the study on the lightest scalar resonances within the 
ENJL model. Our aim is to explore whether one can explain within the model, at 
least at qualitative level, the masses and widths of a(600), J<"o(800), ao(980) and 
/o(980) simultaneously. We find that the model encounters serious difficulties for its 
own reason, though it can be finely tuned, unnaturally, to explain the masses and 
widths of /o(600), J<"q(800) and a (980). However in such a case the ENJL model 
can no longer explain the vector meson spectrum. Furthermore, if we only focus on 
the scalar sector and ignore the problem with vector meson mass spectrum, there 
still remains a problem that the /o(980) is not possible to be described as a member 
of the scalar octet: it has a too small mass. 

This paper is organized as following: section 1 is the introduction, in section 2, 
we make a short review on the ENJL model, especially those materials being used in 
this paper. In section 3, we reconstruct the mass relations of scalar mesons and also 
discuss the tree-level decay widths of scalars. In section 4, using K matrix method we 
construct a unitarized scattering amplitude and find pole locations in each channel, 
numerical results are listed and discussed. We also discuss the N c dependence of pole 
trajectories. Our discussion is also slightly generalized by including the unitarization 
approximation of the more general resonance chiral theory. In section 5, we draw 
our conclusions on the nature of light scalars, based on our study on the K matrix 
unitarization of the ENJL model amplitude. 



2 The ENJL model 

This section reviews how to derive an effective meson chiral lagrangian, involving 
both scalars and vector mesons, from a four fermi interaction. Combining with 't 
Hooft's interaction lagrangian the ENJL model provides the basic tool for our study. 
The method introduced in this section is standard [2, 3, 4]. 



2.1 Bosonization of the ENJL model 

We start from the four quark interactions 

£qcd - ^m-M)q + £^ + J%£ + 0(l/A$, (1) 



with 41 = w 



N r At 



and 4j A l = E KMl) UlI^l) + {L -> R)] , 

c x i,j 

where i,j are flavor indices, ^r,l = (1/2) (1 ± 75) \I> an d the couplings Gs and 
Gy are dimensionless quantities. We adopt the same symbols and definitions as 
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in the third reference of Ref . [ i] . We introduce three complex 3x3 auxiliary field 
matrices M(x), L^(x) and R^(x), which under the chiral group G = SUl(3) x SU r (3) 
transform as 

M -> g R Mg[, 
-> g L L^g ] L and -> g R R^g R . (2) 

By polar decomposition 

M = UH = (3) 
with U unitary, H (and H) hermitian and 

£($) - )^($, ^.fl) = (4) 

where £($)£($) = U. From the transformation laws of M and £, it follows that if 
transforms homogeneously, i.e., 

H^h($,g LtR )Htit(*,g L , R ). (5) 

We can reconstruct the vector fields 

Wt = £L lt ?±?R fl t. (6) 

The transformation properties is 

W£^h@,g)Wfht@,g) • (7) 

After some deduction, one obtains in the Euclidean space the effective action 
V e ff(MQ, £, a, v, a, s,p) in terms of the new auxiliary field variables and in the 
presence of the external field sources v^, a M , s and p [1], 



e 

,2 



r e //(MQ,£,o-,W M ;u,a,s,p) 



exp (- J d x { s^k) trH + uF^Z»WW + w ^ 

J VQVQexp J d 4 xQV E Q, (8) 



where T>e denotes the Euclidean Dirac operator 



V E = 1, V M - i(S - 75 A) - H(x) (9) 



with V M , the covariant derivative 



v M = 0„ + r„ - l 75 (e, - wH) - (io) 



and 

E = gMg + £M*£ (11) 
A = - £M f e . (12) 

The quantities and £ M are those 

IV = \{?[dr - i(v» + a„)]£ + £[d» - i(v, - aj]?} , (13) 

and 

^ = - «'K + Oi£ - eft. - - a M )]e f } = = C • ( 14 ) 

The effective action is in the basis of constituent chiral quark fields Q, 

Ql = iqL , Ql = ql^ ; Qr = SUr » Qr = QrH • ( 15 ) 

With 

exp J VQVQexp J d A xQV E Q = detV E , (16) 
we can get the effective action 

T eff(M Ql f , a, W^; v, a,s,p) 

( N A 2 N A 2 1 
d 4 x < — c x trH 2 H — — tr(W + W + + WW 

+ ln(det|P s |). (17) 

Using proper time regularization and heat-kernel expansion method [18], we get an 
effective Lagrangian of meson fields from the ENJL model. One can also use other 
regularization method [19] to get similar effective lagrangian. 

2.2 Gap Equation and the 't Hooft Interaction 

Here, we are looking for translational invariant solutions which minimize the effective 
action, i.e., 

Jjj \L fl =R IJ ,=0£=l,H=<H>;v fl =a, ll =s=p=0 ~ U ; l i8 J 

where < H >= diag(M u , Ma, M s ). The minimum is reached when all the eigenvalues 
of < H > are equal, i.e., 

< H >= M Q 1 (19) 
and the minimum condition leads to the so called gap equation 

f N C A 2 f 

J d xTl{x\D E \x)\ L ^ =R ^ =Q £ = i tH= M Q ;v^=a^=s=p=Q = ~4Mq ^q^Q^ \ J ^ X > 

(20) 
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where Mq is the constituent quark mass. From Eq. (20) one further gets, 



Gc MA 

M Q = lSn-l^)M 3 Q , (21) 

x x 



where T(—l,x) denotes the incomplete gamma function 

r(n-2,x = -g)= -e-^- 2 ; n = 1, 2, 3, ... . (22) 

A 2 Jm*/ai z 

The Eq. (21) is obtained using proper time regularization method used in this paper. 
The gap equation (21) is obtained without introducing the current quark masses. 
We can introduce the current quark mass through the external source field s = 
diag(m M , md, trig), m u ^ )S is the current quark mass. Unlike the method used in 
Ref. [10] we just use the gap equation without explicit 577/(3) breaking to avoid the 
complicated calculation in heat-kernel expansion. We need to shift the the SU(3) 
singlet field o"o and the octet filed as again in the broken phase to get the physical 
fields with zero vacuum expectation values in the effective lagrangian. 

The next step is to add the 't Hooft interaction [7], (3(detH + detH* 1 ), where j3 
is a constant characterizing the strength of the anomaly contribution. We get the 
modified gap equation, 

5T' eff (H, ...) 6F eff (H,...) 5{detH + detW) 



\H=<H> — FT7 \H=<H> 



+ P- 



=► -^-r(-l,x)M 2 = — ^-2. (23) 

2.3 The Effective Lagrangian and its couplings 

In the ENJL model, we have six input parameters: 

G s , G v , A x , m q , m s , and (3 . (24) 

The gap equation 

A y s 9 8ir 2 (3M , . 

-Z-r(-l,x)MZ = —gr-2 (25) 

introduces a constituent chiral quark mass parameter Mq, and the ratio 

x = (26) 
We can replace the parameters Gs, Gy and A x with x, Mq, and 

9A = l + 4G v xT(0,x) ' (27) 
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characterizing the 7r - A\ mixing. In the limit of Gy —>■ (gA — ► 1), the ENJL 
model goes back to the NJL model. 

The effective Lagrangian can be written down in the form: 

Ceff = \fl [tr {D,UD^) + tr ( x rf + rf x )] 
~ti [V, U V^ -2M V V,V»] 
-itr [A^ -2M\A^\ 

[/vtr (V/f ) +ig v ti + f A tr {A, V D\ 



2V2 

+C S 

+Li (tr D^D^U) 2 + Z 2 tr (D^U ] D u Uti D»U ] D U U) 
+L 3 tr {D ti U^D^UD v U t D u U) 

+L 5 tr [D^D^U { X ] U + C7+x)] + £7 [ tr (£7^ - X^)] * 

+L 8 tr ( X W^ + X^W) + #2 tr ( x t x ) , ( 2 * 



where 



£5 = \tx (d,Sd»S - M 2 S S 2 ) + c m tr (S X + ) + c d tr (S^) + - (29) 

The relevant coupling constants in above lagrangians are listed in the Appendix A. 
We do not integrate out the heavy resonances (vectors, axial-vectors and scalars) to 
get the low energy constants Li, therefore we use Lj to distinguish them from the 
LECs in chiral perturbation theory [20]. 

Attempts have been made in expressing the low energy constants in terms of 
QCD operators [21]. Nevertheless most reliable estimates and determinations at 
this stage are from phenomenological studies [22]. The couplings fy and gy can be 
then determined from the decay p° — > e + e~ and p — > tttt respectively, with the result 

\f v \ = 0.20 and \g v \ = 0.090. 

The decay a± — > 777 fixes the coupling 

\f A \ = 0.097 ±0.022. 

For the scalar couplings c m and q, there exists controversy due to the lowest scalar 
multiplet. One can take the scalar multiplet including ao(980) as the lightest scalar 
nonet as in Ref. [22]. Using the a — > tjtt decay width and assuming the scalar 
saturation of L5 and Lg to determine c m and c<j, in this way one computes with 
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M s = M ao = 983 MeV, 

\c d \ = 32 MeV, 

\c m \ = 42 MeV, 

and CdC m > 0. (30) 

Alternately, the authors of Ref. [23] consider the scalar multiple to be around 1.2 ~ 
1.4 GeV as the lightest scalars in the L s and gives the value q = c m ~ /,r/2. 

The coupling constants, c m , q, fy, /a and GV have been given from the ENJL 
model in Ref. [24]. As will be shown later, the most important parameter appearing 
in this paper is the axial vector coupling g A . The preferred value of g A is found 
to be around 0.6 in Ref. [21]. In Ref. [25], an estimation gives qa = \ under some 
additional theoretical constraints. 



3 Scalar Mass Spectrum and Decays of Scalar Mesons 

3.1 Scalar Mass Spectrum 

The scalar nonet is denoted as the following, 



/ UlL _|_ gQ I CT 8 
V2 ^ V3^V6 



S(x) 



a 



£i. 



K 1 
K° 

I 2 



75 ~ v * as 



\ 
J 



(31) 



Mass relations for the scalar nonet can be extracted from the effective lagrangian. 
Firstly for the charge and flavor neutral scalars there is a mixing term: 



C 



mixing 



M Vo + Mo 8 a a 8 + Mf 8 <7g 



(32) 



One finds 



where 



and 



M, 



oo 



-(2M 2 + M 2 J - g A M t 



M 2 88 = -(AM 2 - M ( 



2 ) 

a )i 



tii 



M 2 08 



m 2 + m 2 , — 2m 2 K 



M(M 2 -M 2 



A(3M t 



+ 



f 2 



M 2 ao ~ 3g A ml + AM 2 + 2 -g A M 2 h , 

2 



M l K ~ 3^ml + 4M^ + -^M^, 



(33) 
(34) 



(35) 
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and Mq is the constituent quark mass. The Eq. (35) is only exact in the leading 
order of cutoff dependence. After diagonalizing Eq. (32) one gets the masses for 
mass eigenstates and the mixing angle, 

M 2 = - 
2 

Ml = - 
2 



2Ml - g A Ml - J(M$ - Mi 8 y + (M 2 8 ) 2 



2M 2 K - g A M t { + ^/(M 2 - Mf 8 ) 2 + (M 2 8 ) 2 
¥(M K 2 -<) 



tan 29 = 9 , 3 v \, . (36) 

From (35) and (36), one gets immediately two sum rules: 

2Ml - Ml - Ml ~ g A Ml, (37) 
Ml -M 2 a ~ 3g A (m 2 K -ml). (38) 

The first sum rule Eq. (37) has been obtained in Ref. [6]. The second sum rule 
Eq. (38) is in qualitative agreement with the results given in Ref. [9]. They are 
the consequence of U A (1) breaking in the ENJL model combined with linear SU(3) 
symmetry breaking terms. If we do not include the anomaly contribution in the 
scalar mass spectrum, just setting M 2 h = in (35), we can find the mixing angle 
9 = (arctan 2v^)/2 = 9 ld ~ 35.26°. Then the a is a pure non-strange state, and the 
a' is purely strange. The scalar masses become 

K = M l = ^ Ami + 4M 2 , 
Ml = 3g A m 2 K + AM 2 , 

Ml, = 2M 2 K - ml . (39) 

Before jumping into more detailed numerical calculations, we can make some 
simple estimates and discussions at qualitative level with Eqs. (37) and (38). The 
first important thing to notice is that, as already emphasized in Ref. [17], the scalar 
masses appeared in Eqs. (37) and (38) are only 'bare' mass parameters appeared in 
the lagrangian, which, when the interaction becomes strong, can be totally different 
from the pole mass positions. The large width of a (or n) is an unambiguous signal 
for a strong airir (or kK7t) interactions. The large widths are quite often ignored 
when discussing the mass spectrum in the literature. It is often attempted to set 
up SU(3) mass relations among pole mass parameters m. However, a light a with 
a mass around 500MeV as a bare parameter appeared in the lagrangian can hardly 
produce a large width, in any model calculations. On the other side, the parameter 
M for a or k can be estimated to be M a ~ 930MeV and M K ~ 1380MeV with sizable 
error bars [17]. Qualitatively speaking the stronger the resonance couples to the im 
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continuum, the larger the deviation is between m and M. Instead of comparing 
different m, one should firstly examine the relations between different "bare" mass 
parameters, M. Since the former quantities associated with large widths are severely 
distorted by the strong couplings to the pseudo-goldstone pairs, it is not suitable to 
use them to discuss the SU(3) mass relations. For example, we have 

m a <m K < m ao , (40) 

but actually the mass relation should be read as 

IGeV ~ M a < M ao < M K . (41) 

The mass of a 1 is unavoidably large, with or without anomaly contributions. If 
we include the contribution of the 't Hooft interaction, taking for example M K = 
1.2GeV, M a ~ IGeV, g A = 0.6, and Mf h = 0.72GW 2 , we get M ao = 1.02GeV, 
M a , = 1.2GeV, 9 ~ 23.59°, and 

/V| / 0.98 -0.20 

\a'J \ 0.20 0.98 J \-cr a J 

from Eqs. (35)-(38). If qa grows larger, the mass of o ! will be heavier. If we neglect 
the contribution of 't Hooft interaction, and taking for example M a = M a ~ IGeV 
and gA = 0.6, we get M K ~ 1.19GeV, M a i ~ 1.35Ge^. If g& grows larger, masses of 
K and a 1 will be pushed higher too. In both cases, with and without 't Hooft term's 
contribution, the a 1 is problematic within the present scenario to be identified with 
the physical /o(980) state, simply because the former is too heavy. As will be seen 
in the discussion given later in this paper, that in order to explain the large width 
of a and k in a dynamical approach, one needs large values of Mq and g A - The 
immediate consequence is that the ENJL model would predict an unacceptably large 
vector meson spectrum and hence fails to give the correct description to the mass 
of p and ai mesons. The reason for this is because in ENJL model the correlation 
between the parameters of scalar sector and the vector seems to be too strong. This 
is not necessary for hadron physics - in resonance chiral theory, for example, there is 
no such strong correlations between the two sectors. Hence we will in the following 
only focus upon the scalar sector and ignore the problem in the vector sector. The 
experience we are going to obtain is still meaningful - if not within the ENJL model 
itself - in a more general background, as in the resonance chiral theory. 

In the scalar sector, as we mentioned above, the problem remains how to identify 
the / (980) resonance. A possible way to solve the a' and the / (980) problem is 
that since the bare state a' is much heavier it may mix with /o(1370), etc. Without 
instanton effects, the a and a' are ideally mixed and the latter is \ss >. When the 
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instanton effects are taken into account, a 1 may contain a sizable \nn > content 
and hence may have a sizable mixing with the heavier scalar like /o(1370), thus 
reducing to some extent its mass. On the other side, one may identify the a' simply 
to the /q(1500) state, since the mass can be quite close to each other and the 
latter is known to be mainly ss state. Then the /o(980) may be considered as a KK 
molecule [27, 28]. Considering the complicated situation about /o(980), a convincing 
explanation to / (980) is out of the range of the present discussion and remains to 
be explored in future. 

There are six parameters (Mq, x, gA, rn q , m s and (3) in the ENJL model under 
investigation, and there are different ways to choose physical parameters to be fit. 
For example, we can fit f n , f K , m n , m K , M th (through M% h = rrv^ + mL — Im^)^ 
m p and the bare mass of a. In the old literature, the bare mass of a (and also ao 
in the absence of anomaly) is typically 500 - 600MeV which is too small. Since 
roughly there is a mass relation in the chiral limit, M a ~ 2Mq [29], it is difficult to 
increase the bare mass of a within the ENJL model. Adding 't Hooft interaction 
term will even further decrease the singlet a mass. So the first thing to be noticed 
is that it is somewhat unnatural to assign a sigma mass of order lGeV in the ENJL 
model. As can be seen from table 1, the x parameter is quite large, which is not 
natural in the cutoff effective lagrangian approach. Another problem is that fy (~ 
0.1) can no longer be fitted well to its experimental value (~ 0.2). Also the current 
strange quark mass gets unnaturally large when gA increases. Furthermore, besides 
these problems, it is clear from table 1 that when gA gets larger the mass of a' is 
enhanced and deviates more and more from the narrow width state / (980). Barring 
this problem, setting M a ~ lGeV, the bare mass of the k resonance is an output 
which turns out to be ~ 1.3GeV and agrees within expectation. Table 1 provides 
several fit values. 

3.2 Decays of Scalar Mesons 

A serious investigation of the scalar mass spectrum unavoidably requires taking 
unitarization into account. But before doing that, in this section we will discuss 
at tree level the decay widths of light scalars, which can be helpful, though very 
rough, in the understanding of strong interaction dynamics behind. For example, if 
the decay width in a given channel in perturbation calculation is small then we can 
judge that the interaction is not strong and the difference between bare mass and 
pole mass is unimportant. If on the other hand the decay width is very large then 
one may claim that the difference between bare mass and pole mass ought to be 
large. In the latter case one has to find more reliable method to handle the strong 
interaction dynamics rather than calculating decay width perturbatively. 
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Table 1: Experimental values and predictions of the ENJL model for the ve 
low energy parameters discussed in the text. All dimensional quantities are in 
except ml h in GeV. 





exp. 


fit 1 


fit 2 


fit 3 


fit 4 


fit 5 




value 












u 


92.4 


92.6 


92.5 


92.3 


92.0 


91.7 


fx 


112.0 


102.0 


106.7 


112.4 


118.9 


136.1 


m n 


137.3 


137.2 


137.2 


137.2 


137.3 


137.3 


m K 


495.7 


495.6 


495.7 


495.6 


495.6 


495.4 


m lh 


0.727 


0.645 


0.44 


0.31 


0.23 


0.16 


Ml 


~ 930 


856 


869 


881 


892 


903.8 


M a 


984.7 


1039 


1025 


1016 


1010 


1004 


Ml 


~ 1400 


1227 


1274 


1330 


1391 


1458 


M a , 


980 


1360 


1456 


1560 


1669 


1784 


X 




0.175 


0.234 


0.295 


0.356 


0.419 


Mq 




397.0 


395.3 


394.1 


393.3 


393.3 


9a 




0.5* 


0.6* 


0.7* 


0.8* 


0.9* 


m q 




4.6 


6.9 


9.6 


12.7 


16.1 


m s 




114. 


172.4 


240.0 


317.1 


403.5 






9.2 


6.7 


5.2 


4.0 


3.1 



(*) Corresponding to bare masses discussed in the text. 
(*) Values of g& are fixed in the fits. 
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Figure 1: The decay widths of a (left) and k (right) as a function of scalar mass M a 
and M K , respectively; for g A =1.0, 0.8 and 0.6. The units are in GeV. 



We use the effective lagrangian to calculate the decay rates of a scalar into two 
pseudoscalars, at tree level. The a decay width is expressed below, 

r = 

1 <T— >7T7r UJ - G— >7TQTro 



H9A(Mi - 2mi) + mi) 2 ( cos 9 + ^=-) 2 , (43) 

where is the scalar meson mixing angel defined by Eq. (36). When is equal to 
8id, the decay width is maximal. From Fig. 1, one realizes that in order to explain 
the large discrepancy between m a and M a , should not be small, for otherwise the 
decay width is small and the interaction will not be strong enough to develop a big 
difference between m a and M a . Especially, from the Fig. 1, we realize that the light 
M a ~ 500MeV and M K ~ 700MeV can not produce large widths. In the SU(3) 
limit, we have g awn = g^Kw and the width is proportional to g\. The only possibility 
in both cases to get a large width is to increase the bare mass parameters. We also 
plot the decay width of a — > Trr} 8 in Fig. 2. The decay width of a is much smaller 
comparing with that of a and k simply because of SU(3) symmetry. See Fig. 2 for 
illustration. Therefore, as revealed by Figs. 1 and 2, the ENJL model does provide 
a possibility in its parameter space to explain the observed scalar spectrum and the 
vastly different widths simultaneously, at least qualitatively. 



4 Pole masses of scalar resonances in the ENJL model 

4.1 The K matrix unitarization and the pole positions 

Certain unitarization approximation is necessary when a large width is involved. 
The unitarization method has been applied to (resonance) chiral perturbation the- 
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Figure 2: The decay width of (do — > ^f]s) as a function of bare scalar mass M ao for 
values of gu=l-0, 0.8 and 0.6 respectively. 

ory amplitudes, and also to linear sigma models in the literature (see for example 
Ref. [31, 23]). To our knowledge, this paper is the first attempt to apply unita- 
rization to ENJL amplitudes. The scattering amplitudes for two pseudoscalars to 
two pseudoscalars are easily obtainable at tree level in the ENJL model, The single 
channel K-matrix unitarization is the following: 

T^tree 

T = — . (44) 

1 - ipT tTee v ; 

We use the K Matrix amplitude determined from ENJL model to search for pole 
positions of scalars, which are not found in the previous literature. The results, 
corresponding to several choices of g\ = 0.6, 0.7 and 0.8, are listed in table 2. In the 
unitarized amplitudes there are actually quite a few poles in each channel, on both 
sheets. Nevertheless in each channel there is only one pole that falls on the real 
axis in the large N c limit which is just the input pole in the lagrangian.* However, 
the results listed in table 2 should not be understood as accurate in any sense. On 
the contrary, it is known that the K matrix results are crude for derivative coupling 
theories [32]. The results given in table 2 only provide a qualitative guide to the 
underlining dynamics: when qa is small the coupling strength between a and nir is 
small and the width of a is also small. The mass of the a found from the unitarized 
amplitude is therefore very close to its input value. However, when increases up 
to, for example, 0.8, the width of a becomes large, and the pole mass m a becomes 
totally different from the input bare mass, M a . 



*Thc er' is very heavy and lies far above the irn clastic unitarity region and hence we do not 
attempt to make any discussion based on the unitarized amplitude. 
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Table 2: Scalar pole masses 





fit 2 = 0.6) 


fit 3 = 0.7) 


fit 4 (g A = 0.8) 


7T7T - 


-> 7T7T 


985 ± 133z 


1161 ± 229z 


544 ± 306z 


ttK - 


-> iiK 


1423 ± 153z 


1629 ± 202z 


801 ± 360z * 


TTT] - 


-> 1X7] 


1030 ± 31i 


1016 ±45i 


1000 ± 60i 



*: The N c trajectory is marginal. 

4.2 Pole trajectories with respect to the variation of N c 

As stated in last section that, all poles listed in table 2 fall on the real axis in the 
large N c limit. t However, there are other poles on the second sheet and it is checked 
that they all go to oo on the complex s plane when N c — > oo. Hence these states 
are dynamically generated. As discussed in section 4.1, for small values of g A (for 
example, fit 2 and fit 3) the a pole from the ENJL lagrangian has a rather small 
width and a large mass around lGeV (as an input), but it is observed that in such 
a case there still exists a light and broad dynamical pole which disappears when 
N c — > oo. This pole, being dynamical, is certainly not the a pole responsible for 
chiral symmetry breaking in the ENJL model, since the latter is well monitored 
and falls on the real axis in the large N c limit. One may even further ask the 
question whether the experimentally observed /o(600) is the a responsible for chiral 
symmetry breaking (In the present situation corresponding to the light and broad 
resonance when g A is large as in fit 4), or a dynamically generated light and broad 
resonance, which is not the a, when g A is small. To understand why there appears 
a 'dynamical pole' we recall that in general the tree level IJ=00 channel hit elastic 
scattering amplitude, in the chiral limit, may be written as 

T -J— 
aA - 16*/? ' 

T 4 = ^yiA , (A = 22L X + 14L 2 + 11L 3 ) 

r 2 3 s 2 9M 4 q 

^^-^ + ^ + ^[1 + ^]), (45) 

where q ~ 0(^/1^) and are obtained by integrating out all resonance fields except 
scalars. The above expressions generally depict resonance chiral theory amplitudes 
including the ENJL model. The pole position in the chiral limit is determined by 

TFor gA = 0.8 the kappa pole trajectory is actually marginal, the k pole will fall on the real 
axis when further increasing gA- 
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the equation 

l + i(T CA + T 4 + T s ) = 0. (46) 

In Eq. (46) if we set T4 and T s vanishing, we get the 'current algebra sigma' pole 
as already discussed in Ref. [32]. The N c dependence of the 'current algebra sigma' 
pole position is y/s po i e ~ \fN~ c an d is ruled out through the study of Ref. [32, 33]. If 
setting T 4 = and \/2c d / f n = 1 in Eq. (45) we recover the linear a model amplitude.* 
In such a case the a resonance is light and broad when the bare mass of o is around 
lGeV. In general, however, if we neglect the logarithm term in Eq. (46), which is 
suppressed when s is large, it is not difficult to show that on the second sheet there 
exists, except the stable pole in the large N c limit, another pole with the property 
s P oie ~ v^Vc on the second sheet of complex s plane. Notice that dynamical pole 
obtained from Eq. (46) contains a different N c behavior comparing with the 'current 
algebra <r': the latter behaves as ^Js po \ e ~ \fW c . The different N c dependence of the 
pole trajectory actually reminds us that the property of the so called 'dynamical' 
pole can be highly (unitarization) model dependent. 

To prove the illegality of the light and broad dynamical pole generated from 
simple K matrix unitarization of the tree level ENJL amplitude when qa is small, 
we make use of the low energy matching method developed in Ref. [34] (see also 
Ref. [33]). For -kit scattering S matrix poles (on the second sheet) obey one relation: 

£sFl<~ 0(iV ' rl) - (47) 

R H 7T 

where M\ and Gr are functions of the pole mass zq of resonance R [14]: 

Mi M = R. W + Im W. Im ';°f°)' , (48) 

Gn(zo) = * m( ? } V] , (49) 
Re[z p(z )\ 

p(z ) = yfl - 4m2/z • (50) 

No matter where on the second sheet does the pole locate, one always has 
Gr/(Mr — 4ml) > 0- The pole solution of Eq. (45) with the property s po i e ~ \/Nc 
corresponds to M R oc O(v^Vc) and Gr oc 0(1). Its contribution to the l.h.s of 
Eq. (47) is 0(l/y/Nc), meanwhile, the right-hand side is 0(1/N C ). The only pos- 
sibility to satisfy Eq. (47) is that the contribution of a such 'dynamical' pole is 
canceled by a spurious pole on the physical sheet, whose contribution is also of or- 
der of 0(1/ y/Nc). This is just the case what we found from solutions of Eq. (45). 



■'■In ENJL model we have approximately V2cd/ fn — Qa- In resonance chiral theory it is found 
that \/2cd/ frr — 0.53 [22] which corresponds to qa — 0.65 here. 
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Hence we demonstrate that the dynamical light and broad pole in the ENJL model, 
generated in the present simple K matrix unitarization, is spurious. 

We can also check the [1,1] Pade amplitude. The pole position of the [1,1] Pade 
amplitude is determined by the equation, in the chiral limit, 

T 2 -(T, + T s -t\T 2 \ 2 ) = 

! A-T^As-2Ml + J^) + ^{-^-f = ^ (51) 



As before we neglect the logarithm term in above. It is straightforward to show 
that if no accidental cancelation occurs, when N c — > oo, there exist two poles on the 
second sheet of complex s plane, one is on the real axis and the other remains on 
the complex s plane: s po i e oc 0(1). At the same time, a spurious pole on the first 
sheet will be found, which is also s po i e oc 0(1). The latter exactly cancels the second 
sheet pole to meet the N c order of the l.h.s of Eq. (47). Hence the dynamical pole 
with Sp i e oc 0(1) found from the amplitude (51) is also a spurious one. Therefore 
the situation as described by Eq. (51) is quite different from the Pade amplitudes 
constructed from pure chiral perturbation theory [36, 37]. The latter is obtained by 
further integrating out the explicit scalar degree of freedom. There one does find 
that the dynamical pole falls on the real axis in the N c — > oo limit. 



5 Discussions and Conclusions 

In this paper we discuss the possibility whether one can understand the light and 
broad o and k, together with the narrow ao(980) and / (980) in a same SU(3) 
nonet, in the ENJL model. We find that the ENJL model is quite reluctant for this 
picture. One difficulty is that the a' resonance is simply too heavy to be identified 
as the /o(980) meson. One has to call for other mechanisms for the rescue. For 
example, the mixing with / (1370) and/or / (1500); or that / (980) is simply a KK 
molecular state [28, 35]. Beside this difficulty, however, the ENJL model can give a 
rough but unified description to the light and broad a, k and the narrow ao(980). For 
sufficiently large qa and an input bare a mass around lGeV, a simple unitarization 
approximation generates a light and broad a resonance. The difference comes from 
the fact that the a couples very strongly to hit continuum, hence its pole location is 
severely distorted. The price paid for this picture is that the g& and Mq parameter 
have to be unnaturally large. As a consequence, the ENJL model is no longer valid 
for describing the vector meson spectrum. However, if we disregard the constraints 
among parameters of ENJL model, the above picture can be realized without any 
foreseeable difficulty in general. 

We also discussed the fate of dynamical poles generated from the simple K matrix 
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unitarized hit scattering amplitude, when g& is small. It was confusing to notice 
that, in such a case, there still exists a light and broad dynamical pole which might 
be identified as the observed /o(600) resonance, besides the input heavy (and nar- 
row) a pole. However, we find that this dynamical pole maintains a wrongful N c 
behavior which has to be canceled by an accompanying first sheet pole, hence vi- 
olating analyticity and should be spurious. The lesson we learn from this study is 
that one has to be extremely cautious when trying to give a physical meaning to a 
dynamically generated resonance pole from a unitarized amplitude. The property 
of the latter can be highly model dependent. Finally further efforts have to be made 
in order to generate successfully the light and broad scalar spectrum from a general 
resonance lagrangian containing both the scalar and the vector sectors. 

Acknowledgement: We would like to thank Zhi-Hui Guo and Juan Jose Sanz- 
Cillero for helpful discussions. This work is supported in part by National Natural 
Science Foundation of China under contract number 10575002, and 10421503. 
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A Effective Couplings in the Effective Lagrangian 



We list in the following parameters of effective meson lagrangian obtained from 
ENJL model. The following expressions are found in agreement with those given 
in Ref. [24]. In the calculation, only those regularization scheme independent terms 
(leading terms in cutoff dependence) are kept. 

16tH w 16tt 2 f 2 T ■ g A 

2 _ 3 A x 1 _a^2 9a 



Mi = - - * „,„ ; = 6M, 

r2 , C if2 



2GV(A 2 )r(0,x) "^ Q l-g A 



Mi = {M^ + 6Mj)/(l 

^ = ^ ^ [(i - <£)r(o, x) + 2^r(i, x 

AT 2 

A ^ = T6^3 [3r(0 ' x) " 2r(1 ' x)] ' 



£ P [r(-i,x)-2r(o,x)] 



Cd = jh^Aim^ -r(i,x)] 

1D7T Ac 



^3 
#2 



1 AL 



2167T 2 

iV c 
16vr 2 

167T 2 



r(o,x 



24 

i - g\f 



-r(o,x) 



24 

:i-^i) 2 " 



16tt 2 
16vr 2 

167T 2 

N c 

167T 2 



[r(o, x)]^m q 

-r(0, x)g A M Q 
T(-l,x 

r(-i,x) 



1 

2B' 

1 



12B 



Mr 



8B 2 
M} 1 



AB 2 



In the following we list expressions of scalar mass and pseudoscalar decay con- 
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stants. Some of them are not found in the previous literature. 



Mt = 



Ml = 



M 



an 



K 



Moo 



Ml 



M 2 08 



Ml 
Ml 



Ml + ^(A(3M Q + -^=(3v - -j=Pvb), 



iv c r[o,Xj 



N c T[0,x 
16tt 2 A2 



4Mj[2 



3 2 1 

12 



Mr 



4MJ + 4M Q (m u + m d + m s ) + -^vqMq 

V3 



/3- 



Vq 



N C 
16tt 2 



T[0, 



.r 



12 

4MJ + 2M Q K + m d + Am s ) + -j=v Q M Q - 2\fcv % M^ 



+ P 



2M Q + + ±v 8 



4:V2(m d + m u - 2m s )M Q 
2v 8 (4V3M Q + 







i 

2 
1 
2 



M 2 + M 
r2 



8 2 s " \A M oo " M| 8 ) 2 + (M 2 8 ) 2 



M 2 + M 8 2 8 + yj (M 2 - Mf 8 ) 2 + (M 2 8 ) 2 
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fl = ^mx]g 1A M" Q + ^T[0,x}g A M Q (^= + ^= + m q ), 

Ik = ^ m q(9ia + 92A) + ^ r [°' x}9aM q (-^=v - -j=v 8 + m s + m q ), 



ml = 2M Q m, 



m K = M Q (m q + m s ) 



v 



''8 



r[-i,s] 
Vr[o,z]' 
r[-i,x] 



9iA = - Trr( 6m 9 + 2Vsv + V6v 8 )(l - g A )), 
92A = g A (l- ^ r (6m s + 2V3v -2V6v 8 )(l- g A )), 

oMq 



2{2m q + m s ) T[-l,x] 
2\/2(m q — m s ) Y[— l,x] 



Uo and i>8 are the vacuum expectation values of the SU(3) singlet field cto and the 
octet field cig of the lowest order of the current quark mass in the broken phase. f n 
and fx will get the tadpole's contribution. 
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